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Abstract 

We study Dirac operator zero-modes on a torus for gauge background with uniform field 
strengths. Under the basic translations of the torus coordinates the wave functions are 
subject to twisted periodic conditions. In a suitable torus coordinates the zero-mode wave 
functions can be related to holomorphic functions of the complex torus coordinates. We 
construct the zero-mode wave functions that satisfy the twisted periodic conditions. The 
chirahty and the degeneracy of the zero-modes are uniquely determined by the gauge back- 
ground and are consistent with the index theorem. 



1 Introduction 



Toric geometry appears in periodic systems. In many interesting cases background gauge 
field is also introduced. We may then consider topologically nontrivial gauge fields. They 
affects the spectrum of the hamiltonian in an interesting way even if the charged particle 
does not touch the magnetic field. In Dirac theory on a compact space we have chiral zero- 
modes for gauge background with nonvanishing topological charge. This is well-known as 
the Atiyah-Singer index theorem [T]. It can be used as a mechanism of keeping fermions 
massless in Kaluza-Klein compactifications. 

If the gauge fields carry topological Information of the periodic system, it cannot be pe- 
riodic but changes by a gauge transformation under translations by a period of the system. 
Similarly the boundary condition for wave functions of a charged particle becomes twisted 
from the naive periodic boundary condition. This makes the task of solving eigenvalue prob- 
lems nontrivial. In this paper we investigate Dirac fields in an abelian gauge background 
of constant field strength on a torus and construct zero-mode solutions of the Dirac opera- 
tor. They are interesting in gauge theories in connection to chiral anomaly |3 E] and also 
in Kaluza-Klein scenario as a mechanism to have massless fermions. In Ref. |lj Sakamoto 
and Tanimura developed Fourier analysis for wave functions satisfying twisted boundary 
conditions on an arbitrary torus. They analyzed the eigenvalue problem of magnetic lapla- 
cian. Their arguments are based on the representation theory of magnetic translation group 
13 E]. Though their approach is powerful and systematic, it is possible to find Dirac 
operator zero-modes in a way that is more accessible to physicists. Namely, we directly solve 
massless Dirac equation on a euclidean torus and find solutions satisfying twisted boundary 
conditions. We show that the zero-mode wave functions can be expressed as gaussian wave 
function times a function of complex torus coordinates. We find that the twisted boundary 
conditions determine the zero-mode wave functions up to an overall normalization. 

This paper is organized as foUows. In the next section we argue general features of gauge 
potentials on a torus. We also show periodicity conditions for wave functions. In Sect. El 
we solve the massless Dirac equation on T^. This system exhibits some characteristics of 
Dirac operator zero-modes on tori. We show that the zero-mode wave functions are related 
to holomorphic or anti-holomprphic functions satisfying periodicity conditions with a twist. 
In Sect. m Dirac operator zero-modes in an arbitrary even dimensions are analyzed for field 
strengths of a Standard form. Most general field strengths are treated in Sect. El The zero- 
mode wave functions can be related to holomorphic functions. We derive the periodicity 
conditions for the holomorphic functions. In Sect. IHl we solve the periodicity conditions 
and give a construction of zero-mode wave functions. Sect. 13 is devoted to summary. A 
periodicity relation for general translations is derived in Appendix ^ In the appendices 
IbI and 10 we give proofs of the key mathematical facts that are used in this paper. The 
normalization constant of the zero-mode wave functions is given in Appendix O 
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2 Gauge Potential on a Torus 

A d dimensional torus T*^ can be regarded as a quotient space R'^/A, where A is some d 
dimensional lattice. Or equivalently, it can be obtained from a d dimensional parallelogram 
by identifying its opposite d — l dimensional faces. In this paper we take a hypercubic regular 
lattice A = Z'^ generated by the Standard basis vectors Ca (a = 1, - ■ ■ ,d) with (ca)^ = 6^ 
and define T'^ = R'^/A'^ by identifying all the points that differ by a vector in A. In Ref. 
^ a more general torus is considered by taking an arbitrary set of vectors in place of the 
orthonormal basis e^. It is straightforward to carry out this kind of generalization in the 
foUowing development. 

Abelian gauge fields on T'^ 



A = ^A„(x)da;" 



(2.1^ 



are topologically classified by a set of integers 0a6 given by 



5a6 



2tt 



(2.2) 



where Dat is an (a, 6)-plane taken in T'^ and F is the field strength 2-form 



(2.3) 



We see that 2TT(j)ab is the flux through Dab- 

That (pab must be an integer can be understood by considering a parallel transport of a 
wave function along a closed curve C taken in Dab as depicted in Fig. [U The phase factor 
arising from the parallel transport along C is given by 



exp 



-t I A 
'c 



(2.4) 



In applying the Stokes theorem there are essentially two different possibilities of choosing 
a two- dimensional surface with C as its boundary. One is the region enclosed by C, the 
shaded region of Fig. ^^a). It is contractible to a point. The other is the shaded region of 
Fig. ^^b). This is also an allowed surface since the left side of the square is identified with 
the right and the upper side with the lower. These two must yield the same phase factor. 
This implies 



exp 



-t I F 



exp[— 2'7r2<; 



1. 



(2.5) 



We thus find that (pab must be an integer. 
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Figure 1: Two possible choices of two- dimensional surface with C as its boundary. The 
square denotes the (a, 6)-plane Dab taken in T*^. 

The gauge potential A cannot be single-valued on T'^ if it belongs to a topologically 
nontrivial sector. However, the field strength F must be well-defined on T'^. In other words, 
it must satisfy the periodic boundary conditions 

Fab{x + ec) = Fab{x). (2.6) 

This imphes that Aa{x) and Aa{x + e?,) can only differ by a gauge transformation as 

Aa{x + eb) = Aaix) + daXb{x), (2.7) 

where A;, may be a function on R'^. This also gives rise to the transformation property of a 
wave function ip under the translations as 

^Pix + Ca) = ga{x)^{x), ga{x) = e'^^^^\ (2.8) 

where we are assuming that the covariant derivative is given by* 

Daij{x) = {da-tAaix))ijix). (2.9) 

In order for the translations ()2.8|) to give a unique wave function under the combined 
translations x ^ x + Ca ^ x + Ca + eb and x ^ x + Cb ^ x + Ca + eb, the ga{x) must satisfy 
the cocycle conditions 

ga{x + eb)gb{x) = gb{x + ea)ga{x). (2.10) 

That these are fulfilled can be seen by considering the parallel transport of the wave function 
along the boundary of the square D as depicted in Fig. |21 The flux through D can be 
computed as 

F= A = \a{x + eb)-\a{x)-\b{x + ea) + \bix). (2.11) 

Jd, JdD 

*We include the coupling constant in the definition of gauge potential. 
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Figure 2: 

In deriving this use has been made of (j2.7j) . Since the flux through D is equal to 27r times 
an integer, we see that the cocycle conditions ()2.10p are satisfied. 

We now turn to gauge potential with uniform field strength. In this case Fab can be 
written as 

Fab = 27i(Pab. (2.12) 
The gauge potential 1-form A in axial gauge is given by 

A = J2 Fabx''dx^ + abdx\ (2.13) 

a<b b 

Since we can freely shift the constant part by 27r by carrying out a gauge transformation on 
T'^ given by g{x) = e^'^*^'', we can restrict to be in the interval O < < 2tt. Furthermore, 
if det Fab 7^ O, we can remove the constant part of Aa by a suitable shift of the coordinates 
X. This is not obvious in axial gauge. We come back to this point soon in connection to 
symmetric gauge. Henceforth, we restrict ourselves to the case det Fab 7^ O and choose the 
coordinates so as for the gauge potential A to be given by 

A = Y,Fabx''dxK (2.14) 

a<b 

This satisfies complete axial gauge conditions 

Ai(xi, • ■ ■ , Xd) = ^2(0, X2, ■ ■ ■ , Xd) = ^3(0, 0,xs, - ■ ■ ,Xd) = ■ ■ ■ = Ad{0, ■■■ ,0,Xd) = 0. 

(2.15) 

The gauge potential (j2.14p satisfies the periodicity 

A{x + ea)=A{x)+dXa{x), (2.16) 
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where Xa{x) is given by 

d 

Xa{x) = J2 F-bx\ (2.17) 
h=a+l 

Due to ()2.12|) the gauge transformation ga{x) in ()2.8|) is periodic in R*^, i. e. , 

9a{x + e},) = ga{x). (2.18) 

We can regard it as a gauge transformation on T'^. The wave function %p must satisfy the 
periodicity ()2.8|1 with Aa given by ()2.17j) . 

Though the axial gauge is useful in our analysis, it lacks the rotational covariance. The 
shortcoming can be remedied by using the gauge potential in symmetric gauge given by 

1 

h=l 

It is related to the potential in axial gauge by a gauge transformation in R*^ 

Aa = Aa + daA, ^ = ^'"^'^-b- (2.20) 

a<b 

The wave function in symmetric gauge is related to the wave function in axial gauge by 

^(x) = e'^^^'Vla;). (2.21) 

The gauge potential ()2.19|) is not of a most general form. We can add an arbitrary constant 
term to the potential without affecting the field strength. It can be removed, however, by a 
suitable shift of the coordinates x if the field strength satisfies det Fab 7^ 0. 



3 Massless Dirac Equation on 

To solve massless Dirac equation in a background gauge field with constant field strength 
we need several machinery. In two dimensions, however, the problem is rather simple but 
its solution is illuminating to understand higher dimensional cases. In this section we give 
a construction of Dirac operator zero-modes on to illustrate some points that also apply 
in higher dimensional cases. 

The Dirac operator in two dimensions is given by 

\d. + idy-i{A^ + iAy) O y ^ ' 
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where a's are the Pauli matrices and the gauge potential for a constant field strength F12 = B 
in axial gauge is given by 



= 0, Ay = Bx. (3.2) 

We can write B = 2ttNv, where z/ is a positive integer and = sgni?. In two dimensions 
the distinction between B and 27rz/ is not essential and the introduction of v and might 
be considered redundant. As we will see in later sections, the distinction becomes crucial in 
the case of general uniform field strength in higher dimensions. 

The zero-mode of the Dirac operator is a solution to the euclidean massless Dirac equation 

Ip'^j = 0. (3.3) 

In terms of components 

- (::) 

the Dirac equation ()3.3p is given by 

{d^ + isdy + sBx)i)s = 0. (3.5) 

The s = ± label the chiral components of tps if we define the chiral matrix by = —io'ia2 = 
(73. We see from (13. 5|) that ips can be written as 

ijs{x,y) = e-^^«^V.(x + M, (3.6) 

where fs is a function of x + isy to be determined shortly. 

Under the translations x — > a; + (a = 1, 2), the gauge potential (j3.2j) is changed by a 
gauge transformation Xa{x) = By5ai- This gives the periodicity conditions on ip as 

^(x + l,i/) = e'^J^7/'(x,?/), iP{x,y + l)=ilj{x,y) (3.7) 

These can be transcribed into the conditions for fs{x + isy) as 

/^(3; + sN + isy) = e2-(-W|^^) /^(a; + ^sy), 

fs{x + ts{y + N)) = f,{x + tsy). (3.8) 

Since fs{x + isy) are periodic under the translation y y + N, we can Fourier expand them 



as 



/.(a: + z.y) = 5^<e2-"(^+^^^). (3.9) 
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The Fourier coefficients can be determined from the first condition of ()3.8|) . They must 
satisfy the recursion relation 

< = e--^^(2"-^)<_,. (3.10) 
By putting n = uk + r with A; G Z and r = 0, 1, ■ ■ ■ , z/ — 1, these can be solved as 

where are arbitrary constants. The fs{x + isy) is now given by 

v—l +00 

fs{x + lSy) = ^clY^ ^-nsN{uk^+2rk)+2n(uk+r){x+isy)_ ^^_^2) 
r=0 k=—oo 

Since B = 27ti'N, the sum defining /_ does not converge if i? > 0. This imphes that c~ 
and, hence, /_ must vanish. There are only z/ independent /+ corresponding to an arbitrary 
choice of c^. We thus obtain u independent zero-modes with positive chirahty. Conversely, 
/+ does not exists for B < O and there are z/ independent giving zero-modes with negative 
chirahty. 

The absence of zero-mode of Di — iD2 for 5 > O can be shown without solving the Dirac 
equation. Using 

{Di + iD2){Di-iD2)=Dl + Dl-B, (3.13) 

we get 

/ d^x\iDi - tD2)^P-\^ = [ d^xi\Di^^\^ + \D2^P-\^ + B\ij_\^). (3.14) 

The rhs of this expression cannot vanish if ip_ ^ 0. This imphes that Dx — iD2 has no 
nontrivial zero-modes for i? > 0. Conversely, Di + iD2 has no zero-mode if i? < 0. 
Combining (j3.6|) and (j3.1!2j) . we finally obtain the normahzed zero-modes 

iprR = u+ipr,+ ix, y), {B > 0) 

^,L = u_^,,_(x,i/), (5<0) (3.15) 
where we have introduced the basis of two-component chiral spinors u± given by 

= • = (i) ■ 

The ipr,s{x,y) are defined by 

-^[u(sNx-k)-r> +2iTis{uk+r)y 



i^rA^jV) = (2z^)^ ^ e 



k=—oo 

(2z/)i Y, e-^(^'^^"-")'+''^*^"^ (3.17) 

n=r (mod u) 



where the last sum is taken over integers of the form n = uk + r {k ^ Z). The overall factor 
is so chosen that the wave functions satisfy the orthonormahty 



d Zlp* ^Ipr' ,s = Sr,r' . (3.18) 

The iprR and tpri are eigenstates of the chiral matrix T^. If we denote the number of 
positive and negative chirahty zero-modes by ria and rii and define the index of Ip by 

index^ = ur — n^, (3.19) 

we find index^ = uN = B/27r irrespective of the sign of B. This is completely consistent 
with the index theorem. 

B \ f 

index^ = — = — F. (3.20) 

2% 2% Jrp2 

4 Extension to Higher Dimensions 

The analysis given in the previous section can be easily extended to arbitrary even dimensions 
d = 2n in the case that the field strength 2-form takes a special form 



F = J2 Bkdx^''-^ A dx2^ (4.1) 



k=l 



where Bk is a flux density through {2k — 1, 2A;)-plane and z/^ = \Bk\/2Ti is a positive integer. 
Most general uniform field strength can be made into the form 1)4.11) by a suitable coordinate 
rotation. In the rotated coordinates, however, the periodicity conditions ()2.8p become com- 
plicated. We postpone the analysis of the Dirac operator zero-modes for the most general 
uniform field strength to later sections. In this section we argue the Dirac operator zero- 
modes for the special gauge background given by ()4.H) and develop some machinery needed 
in later sections. 

The 7-matrices {(^ = 1? ■ ■ ■ id) 'm d = 2n dimensions can be obtained recursively by 
the relations 

ri =ai®l2n-i, =a=^®l2„-i, r^ = a3®r:i2, (a = 3,--- ,c/), (4.2) 
where Iat is the N x N unit matrix and F" are 2" x 2" matrices. Explicitly they are given 

by 

p2fc+l ^ ^3 (g) . . . ^ ®l2 O • • ■ O lo, 



'-2 

-v ' ^ ^ 



k times n-fc-l times 

^2k+2 _ „3 ^ ^ „3^ ^2 , 

-v ' < ^ 

k times n-k-1 times 



Tl^+^ = a^0...0 (T^ ®l2 ® ■ ■ ■ ® I2 . (4.3) 
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The chiral matrix F^+i is given by 

r,+i = H)" ^^ ■ ■ r^^ = a'^-y^al (4.4) 

n timcs 

The gauge potential that gives rise to the field strength ()4.1|) is a simple extension of 
(13. 2p and is given by 

A2k-i = O, = B,x'"'-\ (4.5) 

The Dirac operator is then given by 

n-l 

IP = T'^ida - tAa) = ^ (T^ ® ■ ■ ■ ® ® {(r^D2k-i + cr^D2k) ® I2 ® ■ ■ ■ ® 12- (4.6) 

The components of the Dirac spinor can be defined by introducing a basis of the represen- 
tation ()4.3|1 defined as direct products of two-component spinors Ug (s = ±) given by ()3.1(i|l . 
They satisfy 

a'^Us = sUs, a^u_ = u^, a^u^ = m_ (4.7) 
with cr^ = |((T^ ± ia"^). The basis of spinors in d = 2n dimensions is then given by 

Usi®---®M«„- (4.8) 
Then an arbitrary spinor ip can be expanded as 

^(^) = X^ «si ® ■ ■ ■ ® Ms„^si-s„(x), (4.9) 
{s} 

where V'si-'-sn stand for the components of ip{x). Using a^Ug = U-s and cr^Ms = isu^s, we can 
write the Dirac equation in component form as 

n 

^ Si ■ ■ ■ Sk-l{D2k-l - iSkD2k)lpsr-s,_,~s^s^+,-sAx) = 0. (4.10) 
fc=l 

These are the generahzation of (|3.5p to d = 2n dimensions. Since each equation contains n 
different components of ip{x), one might think it difficult to solve them. As we will see soon, 
this is not the case. We show that each component tpsi---sn must satisfy uncoupled first order 
differential equations similar to (j3.5p . 

To show this we first compute the square of the Dirac operator. By noting the commu- 
tation relations 

[D2k-u D21-1] = [D2k, D21] = o, [D2k-i, D21] = -iBJku (A;, / = 1, ■ ■ ■ , n) (4.11) 
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we obtain 

n 



k=l 
12 _ n2 



where we have introduced the Laplacian D = . In the present representation of F" 

we have 



p2fc-ip2fc = ^ . . . ^ I2® ®l2 ® ■ ■ ■ ® I2 . (4.13) 



v— 

k—l times n— k times 



This gives 

n 

ip^i) = ^Us,<^---<^Us„(p^ + ^SkB^i)s^...s^. (4.14) 

{s} k=l 

We thus find that each component of the Dirac operator zero-modes must satisfy 

n 

(^D^ + J2^kBk)iJs,-s„ = 0. (4.15) 

k=l 

By multiplying to this expression and then integrating over T'^, we obtain 

/ d''xy^\{D2k-i+iSkD2k)i^s,-sS = 0, (4.16) 
k 

where use has been made of the relation 

{D2k-i - iskD2k){D2k-i + iskD2k) = Dlf^_^ + DIj^ + SkBk- (4.17) 

We thus find that V^si-'-sn must satisfy 

{D2k-i +iSkD2k)^ s,-s„ = 0. (4.18) 

These are essentially the same with the two- dimensional equations ()3.5|) . As we have shown 
in the previous section, D2k~i + iskD2k has no zero-modes if sgnskBk < 0. This imphes that 
'ipsi-sn = O unless Si = sgn_Bi, ■ ■ ■ , s„ = sgni?„. The zero-modes of Ip is then given by a 
tensor product as 

ij = iJl{Xi,X2) <S)1p2{x3,X4) ® ■ ■ ■ ^iJn{Xd-l,Xd), (4.19) 

where ipk is a zero-mode of a^D2k-i + o'^D2k- There are z/^ zero-modes of a^D2k~i + cr'^D2k 
with the two-dimensional chirahty sgni^^. The chirahty of ip is the product of each chirahty 
of the two-dimensional spinor ipk and coincides with sgnBi--- Bn- Since the number of 
independent zero-modes of Ip is given by z/iZ/2 ■ ■ ■ t'n, we obtain the index relation 
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5 Extension to Arbitrary Field Strengths 



In the previous section we have considered rather special gauge field that can be treated by 
the method developed for the two dimensional system described in Sect. El A general uniform 
field strength can be more complicated. The fiux through (z,j)-plane may be nonvanishing 
for any i ^ j. Fortunately, an arbitrary uniform field strength Fab can be made into a 
Standard form by an SO((i) rotation as 

/ O 

-5i 



c,d 



cdi 



B = {3, 



ab) 



O 



O 

-B2 



B2 

o 



v 



o 

-B„ 



Bn 
O 



(5.1: 



where R = (Ca) is a real orthogonal matrix satisfying RR^ = R^R = 1 and det R = 1. Since 
we are assuming det Fab 7^ O, any eigenvalue B^ is nonvanishing. This ensures that the Dirac 
operator has a nonvanishing index. The essential difference between (|4.H) and (jS.ip is that 
the eigenvalue Bk in the latter is not necessarily a 27r multiple of an integer. We also note 
that the field strengths Fab can be expressed as 



ab 



/■2k-l/-2k 
^b 



cf). 



(5.2) 



We will frequently use the relation in the foUowing development. 

The goal of this paper is to find the solutions to the Dirac equation coupled to the gauge 
potential in axial gauge ()2.14|) . It is given by 



a b<a 



What makes the task nontrivial is the periodicity conditions 
are given by 



1p{x + Ca) 



(5.3) 

In the present case they 
(5.4) 



In Ref. jl] these are referred to as twisted periodicity conditions. We will see that they 
essentially determine the solution. 

In the axial gauge, however, the rotational symmetry is not manifest. So, we carry out 
the gauge transformation ()2.20j) and move in symmetric gauge. The Dirac equation ()5.3|) 
then becomes 



$^r«(5a + ^5^i^a6X^)V;(x) = 0. 



(5.5) 



12 



We now introduce new orthogonal coordinates y" by 

?/'^ = E^'C (5.6) 

b 

The Dirac equation (|5.5|) can be converted to the form 

a ^ b 

where '^{y) = Sip{x) with S'^T'^S = ^^T^i^^. The S stands for the spinor representation of 
the SO((i) rotation R. 

We again go back to axial gauge in the coordinate system y. This can be achieved by 
the inverse of the gauge transformation ()2.2H) . The Dirac equation ()5.3|) finally becomes 

d 



j:^idZ-a+'j:Baby'')ny) = o, (5i 

^ b<a 



where \l/ is related to the original ip by 



The Dirac operator on the rhs of (j5.8|) is of the from given by ()4.(i|l . This imphes that the 
Dirac equation can be reduced to the form 

+ ^^^^ + ^kB,y''''\ ^^1-. =0. (fc = 1, ■ ■ ■ , n) (5.10) 
We see that can be expressed as 

^s,-sM = e"^^^^''^^(^"")V....s„(^), (5.11) 

where fsi --s„ depends only on the n complex variables z'' = y^fc-i _|_ isj^y'^'' (k = 1, ■ ■ ■ ,n). 

We now turn to the periodicity conditions for the wave function The translation 
X ^ X + Ca corresponds to a translation of ?/ by a vector (a = (Ca) '^^^ ^e seen from 

y' = E^'C' - y" = + ^^yC = {y + Caf. (5.12) 

c c 

Using the periodicity ()5.4|) under the translation a; — > x + e^, we can find the periodicity of 
under the translation ()5.12j) as 

■^{y + Ca) = e*S''<=^'^''^(^'+^^")^(?/). (5.13) 
The overall factor of gauge transformation 
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satisfies the cocycle conditions analogous to ()2.10|) . This can be shown explicitly by noting 
the relation ()5.2j) . The ga{y), however, is not single-valued on T'^. This is contrasted with the 
corresponding factor appearing in ()5.4|) . where the gauge transformation is periodic under 

Under the translation ()5.12|) the are transformed as 

,>^^,'' + E'^ with E'^ = C-' + ts,C. (5.15) 

The periodicity conditions ()5.13|) can be transcribed as the conditions for fzj^...s^{z). They 
are given by 

fs.-sA^ + ^a) = eS^-^^^^^"^(^^+^^^)^...„(^). (5.16) 

These ensure that the wave functions defined in R'* can be regarded as those on T'^. We 
will see that they imposes very tight constraints on the wave functions and are sufficient to 
determine fs^-s„- 

Before closing this section we generahze ()5.1(ij) to an arbitrary translation by a vector 
= Yla-^"''^a,, where A^" (a = 1, ■ ■ ■ ,2n) are integers. This can be achieved by applying 
()5.16|) repeatedly to fs^...sSz + S) until we arrive at /^^...^^(z). Since we will use a similar 
formula in later sections we give the periodicity conditions for an arbitrary translation in 
Appendix^ They are given by 

^...„(^ + H)=e^^-^^^"^^^(^^+H-fi:.<.^"^^^-/,^...,„(.), (5.17) 

where ^^fc-i ^ ^la^^C^"^ ^^^^ P^^^ °f '^a- 

The periodicity under the translation 2; — > z + S becomes much simpler if we can find a 
set of integers N"" for which C^'^"^ {k = 1, ■ ■ ■ vanish. In the next section we show that 
we can find a set of translations under which fsi---s„{z) is periodic or anti-periodic. 



6 Solving the Periodicity Conditions 



In this section we give a construction of /si...s„(;z). To achieve this we search for a basis 
vectors other than {ca} for the lattice Z'^^ and convert the periodicity conditions (j5.16p into 
a tractable form. 

We first note that the field strengths Fab can be made into a block- diagonal form 

/ O z/l \ 

-ui O 

O z/2 

-1^2 O , (6.1) 



1 



l^ab — T^-^a J^b ah 
ZTT 



O Vn 
-Vn O / 
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where L = {L J') is a regular 2n x 2n matrix with integral elements satisfying | detL| = 1 
and z/p (p = 1, ■ ■ ■ , n) are positive integers with the property that Up divides Ug for p < q. 
As before, we are assuming a general field strengths Fat with det Fab 7^ O and any z/p is 
nonvanishing. For given Fab we can find L in a systematic way. To illustrate this we give a 
proof of this lemma jH] in Appendix iBl 

We denote L by 2n vectors Mp = {MP'") and Np = {Np") as 

L = {M,,N,,--- ,M^,N^). (6.2) 

Then the relation can be expressed as 

F(Mp, Ng) = Up5pg, F{Mp, Mg) = F{Np, Ng) = O, (6.3) 

where we have introduced a skew symmetric bilinear form 

^(^'^) = ^E^'^''^V- (6.4) 

a,b 

It is integer-valued if both ^ and rj are integral vectors. 

As we show in Appendix ^ any vector in A = Z*^ can be uniquely expressed by a linear 
combination of the 2n vectors Mp and Np with integral coefficients. This implies that A can 
also be generated by the basis {Mp,Np}. In defining T^" = R^"/Z^", the 2n dimensional 
hypercubic region O < a;" < 1 in R^" can be replaced with the rectangular region D given 

by 

x^ = ^(Mp'^e'' + ApV)- (O < e, < 1) (6.5) 

p 

Using ()6.3|) . we can solve this with respect to and 77^ as 

^ Fix,Np) ^ ^ ^ F{Mp,x) _ 



The volume of D can be easily found as 



[ d^''x= I I c^"r7|det(Mp^A^p'^)| = 1, 

J D Jo Jo 



(6.7) 



where use has been made of | det L\ = 1. These imply that and 77^ can be regarded as the 
coordinates of the 2n-torus T^". Later in this section we consider the translations 

z" ^ z" + Mp"", Np". (6.8) 

These correspond to unit translations of the {C,, 'Z7)-coordinates. 
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In Sect. EJ we have used the fact that the field strength can be made into a Standard 
form ()5.H1 by an S0(2n) ration R = ((^^). The choice of R, however, is not unique but we 
can replace Cl^~^ and Ca^^ by 



C' cos 4 + c sin e,, -C-' sin 9, + C cos 6,, (6.9) 

where 6k is an arbitrary rotation angle. By utilizing this arbitrariness we can always choose 
to satisfy 

E^Af-' = o, (A;,p = l,--- ,n) (6.10) 

a 

where Np = {Np"-) are those vectors that constitute the L as given by ()6.2j) . In Appendix \C\ 
we give a proof of the existence of such a basis C* = {(^). 

For the choice of satisfying ()6.1()|1 the periodicity ()5.16|1 under the translations 

z z + np with fij = isfc E A/'/Ca'^ (6.11) 



reduces to 



where e» is defined by 



U...,„{z + np)=e-''^^-U..s„{z), (6.12) 



^-Y^Np'^Np'F,,. (6.13) 



There are exactly n such translations corresponding Np {p = 1, - ■ ■ ,n). Since is an integer, 
we see that /^^...^„(z) is either periodic or anti-periodic under the translations (16.111) . 
It is convenient to introduce new variables by 



z' 



5] 5^ .,iV/C>^. (6.14) 



p a 



By noting ()5.2|) and ()6.10|) . this can be made into the form 

J2skBkC-'z>' = 2nJ2sapW'', (6.15) 



k 

where Sap is an integer defined by 



Sap = 7^J2^-''NpK (6.16) 



2lT 

b 



From ()6.3|1 we see that Sap satisfies the orthogonality relation 

Y,Mp''sga = i^pSpg. (6.17) 
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This enables us to solve ()6.15|) with respect to as 

^Y.Y.s,B,M,'^C-^z\ (6.18) 



27rz/,„ 

f k a 



The translation z ^ z + flp corresponds to a unit imaginary translation in w 

w ^ w + icp, (6.19) 
where Cp = {6'^) is a unit vector. This can be seen easily as 

z'k = J2J2-^kN,Xl\w'' + t6^p) = z'' + n';. (6.20) 

q a 

If we denote /^^...^„(-z) by (^^^...^„(w), the periodicity (j6.12j) can be written as 

gsr-s„{w + icp) = e-'^'^gs,...sAw)- (6.21) 

Since the change of variables (j6.14j) is nonsingular, the shift oi w ^ w + Ta under the 
translation z z + is uniquely determined as 

= ^ E E SkBkMp\l'-'Zl (6.22) 
P b k 

We thus obtain from ()5.16j] the periodicity relations of gs^-s,^{w) under w w + 'Ta 

gs,...s,Xw + TJ = e^-^^^-(-''+^^") (u;), (6.23) 

where use has been made of fl6.15|) . 

We can generahze (|6.23p for arbitrary T = M^T^ as in (j5.17|) . It is given by 

{w + T)= e'- S. ^^"-^ (-^+^^^) +i '''''''^^gs....s„ (w), (6.24) 
where M" are arbitrary integers. In particular under the translation 

w^w + pp with pp" = Mp^Tl, (6.25) 

a 

where Mp = {Mp"-) is a vector appearing in L, the wave function is transformed as 

9s,-sSw + Pr>) = e2'^^''(-''+^^'''')+^-'^^^^,,...,„(«;), (6.26) 
In deriving this use has been made of ()6.17|) . The ip is defined by 

^P = ^Y1 Mp'^Mp'Fab. (6.27) 

P a<b 

Note that only Wp appears in the exponent on the rhs of ()6.26p under the translation ()6.25|) . 
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We are now in a position to solve the periodicity conditions ()6.21|) and ()6.26|) . From the 
periodicity ()(i.21|l we can express 5'si-- s„(w) in the Fourier expansion as 

9s, -sA^) = E c;-■-e2-S.('^^-^^^)-^ (6.28) 

{rip} 

where {up} stands for the abbreviation of ni ■ ■ ■ n.„. 

The condition ()fj.2f)j) can be translated into the recursion relations for the Fourier coeffi- 
cients. They are given by 

Si---s„ _ si---s„ g-TT^]^ pp''(2n,-e,)+7rt'pPpP+i7rfpi/p (6 29) 

{nq} {nq-UpSpq} ' V'/ 

It is straightforward to obtain the general expression for C|^"'|". We find that c'^^„"y" for 
Ug = Ugkg + V g wlth tg E Z 9.10.6. T g = , 1 , " " " ,h'g — 1 is glveu by 

^SfSn _ ^si---s„^~TTYlp,jPp''iygkpkg-nY^p^^ppikp(2rg-eg)+iTTY^pepUpkp (6 30) 

where c^^!'',." are arbitrary constants to be determined later. In deriving this we have used 
the relation 

PpV, - p/^^p = -iF{Mp, Mg) = 0. (6.31) 

To each component of the wave function labeled by {s k} we have z/i ■ ■ ■ z/„ independent 
Solutions to the recursion relations ()6.30p . They do not, however, necessarily correspond to 
nontrivial zero-mode solutions. The Fourier expansion ()6.28|) must converge. To establish 
the convergence of the rhs of ()6.28|1 . we must show that c^j'.'.'.'^"^ — O as |n| — oo, where 
\n\ = ^yn^~^^^^—^^~n^. It suffices to verify the convergence of the most dominant part 

|e-'^Ep.,Pp''i'gfcpfcg|2 ^ g-Efc^fcSfc(EpE„Ca'''Mp«fcp)'_ |^g_22) 

Since we have det J2aCa^~^^p"' 7^ O' does not converge as \n\ — > oo unless SkB^ > O 

for k = ,n. This implies that the coefficients c^l'.'.'.r" must vanish if {si, ■ ■ ■ ,s„} 7^ 

{sgnSi, • • • ,sgni?„}. We thus find that the only nontrivial components of the zero-modes 
is \l/si...s„ with Sk = sgn(_Bjt). Hence the number of independent zero-modes is equal to the 
number of arbitrary constants c^^'.'.'f^.'^ and is given by 



^'^^i 7^-^"M -/o \« ,\(^aibva„b„Faibi- ■ ■ Fa„b„\, (6.33) 



^ \^27r (27r)"n! 



where ea^bi-- a„b„ is the Levi-Civita symbol in d = 2n dimensions. Furthermore, the chirality 
of these zero-modes coincides with si ■ ■ ■ s„. These are perfectly consistent with the index 
theorem dOO)). 
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We now turn to the wave function ()5.1H) . To write it in terms of the present coordinates 
let us denote w'^ by its real and imaginary parts as 

= uP + ivP. (p = 1, ■ ■ ■ , n) (6.34) 

Then we have from (j6.18|) 

^ 2nJ29pguPu'^, (6.35) 

k p,q 

where gpg is given by 

a,b k 

It reminds us of a metric in geometry. To find the inverse of gpg we note the relations 

E (E ^'^^/Cf ) (E SkBkMX'-'^ = 2vrz//,„ (6.37) 

k a b 



p P a 



These can be verified by using ()5.2j) . ()6.3|1 and ()6.1()j) . The inverse g^"^ can be obtained from 
(IF^F)]) and (OH|) as 

= EE^Cf-^Cf-^M/M,^ (6.39) 

a,b k P'? 

We see from (I6.22j) and (j6.25p that g'^'^ is the real part of pp'^/up. We write as 

Pp' = Mg^' + ^lP'), (6.40) 

where 7^"? is defined by 

Y' = E E BkiC-'d' + Cf -'Cf )M/M,^ (6.41) 

^^^P^l a,b k 

The wave function ()5.11|1 for a given set of {r^} can be written as 

(6.42) 

where the sum with respect to rlp (p = 1, n) should be taken over the numbers of the 
form fip = Upkp + — |ep {kp = O, ±1, ±2, ■ ■ ■ ). We have suppressed overall normalization 
constant in ()6.42|) . 



19 



The wave function ()6.42j) is not the final result. We must go back to the wave function 
ip{x) in the original coordinates. As we have mentioned at the beginning of this section, 

and T]^ defined by ()6.5|) are natural coordinates on T^"'. The coordinates u'^ and are 
related to and rj^ by 



J2 9'%^', vP = r^P + J2 7"V,e. (6.43) 



The translations ()6.19|) and ()6.25|1 correspond to 

^<?_^9^ + under w'' ^ w" + i6^p, 

^"^^" + 51, 7]«-^//« under u;^ ^ w« + p/. (6.44) 

These can be understood by recalling the fact that the translation ()6.19|) and ()6.25p . respec- 
tively, correspond to + Np^ and x" — * a;" + Mp". We thus find that (j6.19|) and (j6.25|l 

are the unit translations of the coordinates and rj^. These are just the basic moves on 
T^". We can also directly show ()6.44|) by making use of ()6.43|) . 

It is now possible to write down the zero-mode wave function ip{x) in terms of and rj^. 
Since the spinor part of \l/ is given by m^-^ ■ ■ ■ Ms„, the normalized zero-mode for a given 
{vp} can be obtained from ()5.9j) and ()6.42|1 as 



o- 

Vdet 



S^Us, ® • ■ ■ ® Ms„Fn...r„(a;)etS''<''^'"''""" , (6.45) 



where F^^'''^^{x) is given by 

f:-:::(x) = e-^s-^^-^''^^M/;-;-(y) 

ni,---n„ 

(6.46) 

In deriving this, use has been made of the relation 

a<b p \ q / 

The orthogonality of ip^^'"^'' with different set of {r^} is obvious by the Standard integral 

f' rf-^e2-^S.^-^''e-2-^^.^^^' = 5^,,^; ■ ■ ■ (6.48) 
Jo 

The computation of the normalization constant is given in Appendix O 

The zero-mode wave function ()6.45p is the generalization of ()3.15|) to arbitrary even 
dimensions. It is interesting to see how the periodicity conditions ()5.4|) . or their generalization 
for an arbitrary integral vector L = L°-ea 

^{x + L) = e'^-<^^^'^''''''ilj{x), (6.49) 
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are satisfied. These can be verified by applying ()5.4p successively to the Ihs of ()6.49|) . Since 
the 2n vectors Mp = {Mp"") and Np = (A'p") span the entire lattice Z^", it suffices to show 
fl6.49|) for L = Mp, Np. In the form given by ()6.45|) it is completely straightforward to check 
these. For instance, under the translation x — > x + Mp the overall exponential factor of il!{x) 
develops an extra phase factor as 

On the other hand is transformed into 

F,7:.::(x + Mp) = e'--P'?''+^-^-p-PF;^i;;-"(a;). (6.51) 

Since ipUp is an integer, the extra phase factors in the rhs' of ()(i.5()j) and (jfj.Slj) cancel each 
other. Similar thing also happens for the translation x — x + Np. 

We finally note that the transformation properties of Fr-^^...r„{x) under ()6.8p become sym- 
metric as 

F:--:"(x + MJ = e'^^'^<bFab{M,^^''-M,''x^+M,^M,'')pr,...r^^^^ 

F:i:::-{x + Np) = e'2T:.<bFM^-''-^.''-^-^.^^p'')F:i:::::{x). (6.52) 

In deriving these we have used ()6.5p . ()6.13|) and ()6.27|) . 



7 Summary 

We have investigated the Dirac operator zero-modes on a torus for a gauge background 
with constant field strength. The gauge potential possesses a characteristic transformation 
properties under the translations corresponding to the periods of the torus. This gives rise 
to nontrivial periodicity conditions for the wave function. In the coordinates where the 
field strength take a Standard biock-diagonal form the Dirac equation for the zero-modes 
can be reduced to first order differential equations for each spinor component. It give a 
kind of holomorphicity conditions and allows us to express the wave function by a set of 
holomorphic functions of the torus coordinates. We have shown that it is always possible to 
find a set of complex torus coordinates in which the holomorphic part of the wave function 
becomes either periodic or anti-periodic under unit translations along the imaginary axes. 
The periodicity conditions determine up to an overall normalization the holomorphic part 
of the wave function. For gauge backgrounds with uniform field strength only left-handed or 
right-handed chiral zero-modes are realized in a manner consistent with the index theorem. 

Throughout the paper we have assumed det Fab 7^ 0. If det Fab = O, then the index 
of the Dirac operator vanishes. We expect two possibilities. One is the occurrence of an 
equal number of positive and negative chirality zero-modes. Another is the absence of chiral 
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zero-modes. As noted in Sect. |2l we cannot remove some of the constant components of 
the gauge potential by shifting the coordinate origin if det Fab = 0. In general the presence 
of a constant gauge potential affects the spectrum of the Dirac operator and eliminates the 
chiral zero-modes. A complete analysis of the eigenvalue problem of the Dirac operator is 
interesting. We will argue this in a separate publication. 

This work is supported in part by the Grant-in Aid for Scientific Research form the 
Ministry of Education, Culture, Sports, Science and Technology (No. 13135203). 



A General Translations 

We first consider a translation z ^ z + riEa with an arbitrary integer n. By a successive 
application of ()5.16|) . we find 

We now consider the most general translation by 

E' = Y,N''Sl (A.2) 

a 

By applying ()A.1|) to 

b 6-1 

a=l a=l 

we obtain a recursive relation 

a=l a=l 

(A.4) 

This gives 

(A.5) 

where use has been made of C^^~^ = ^la^^Ca'^^^- By noting the relation 

a a>b a<b 

(A.6) 

and using ()5.2p . we see that ()A.5|) can be written in the form ()5.17|) . 
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B Proof of the Lemma 



In this appendix we give a proof of the basic lemma mentioned in Sect. El It can be stated 
in the following form; 

Lemma: Any nondegenerate antisymmetric matrix (f) with (pab = ~4'ba ^ Z can be converted 
to a block- diagonal form 

/ O -z/l \ 
z/l O 

O -z/2 

z/2 o 



La'^Lb'^cpcd = ^ab, with z/ 



v 



o 

z/r, 



-z/„ 



(B.l) 



where L is a 2n x 2n matrix with L J' G Z and | det L\ = 1. The z/i, ■ 
of positive integers with z/p+i/z/p e Z. 

To prove this let us consider a space Z^" spanned by the vectors 

e = $^rea, r e Z). 

a 

Then naturally defines an integer-valued skew symmetric bilinear form 



z/„ are a sequence 



(B.2) 



(B.3) 



Let Ml and A''i be such a pair of vectors that z/i = 0(Mi, A''i) becomes the smallest positive 
integer. Then we see that 0(Mi,^) and (j){Ni,^) for any integral vector C, must be divided 
by z/l. This can be shown as follows: Let ^ be such an integral vector with 0(Mi,^)/z/i ^ Z. 
Then we have 0(Mi, kNi + ^) = A;z/i + 0(Mi,^) for any integer k. Since |0(Mi,,^)| > z/i by 
assumption, we can always find a suitable k such that O < 0(Mi, A;A''i + ^) < z/i is satisfied. 
This contradicts the assumption that ui is the smallest positive integer taken by 0. Hence 
0(Mi,^) must be divided by z/i. This also applies to (f){Ni,C,). Using these properties, we 
can decompose any integral vector ^ in the form 



where ^' satisfies 



0(Mi,O = 0(r,iVi 



0. 



[B.4) 



(B.5) 



We next consider the space of integral vectors ^' satisfying ()B.5|1 and evaluate the bilinear 
form 0. Let M2 and N2 be such a pair of vectors that z/2 = 0(M2, A'"2) becomes the smallest 
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positive integer on this restricted space. By assumption we have ui < z/2. To see that 
z/2 /z/l G Z we consider 

(j)(kMi + M2, A^i + N2) = kui + z/2, (B.6) 

where k may be an arbitrary integer. Since kvi + z/2 cannot be a positive integer smaller 
than z/l for any integer fc, we see that z/i must divide z/2. We can also show that 0(M2,i^) 
and 0(A^2,O fo^' any vector S, G Z^" can be divided by z/2. 

This procedure can be continued until we arrive at M„ and Nn with = F{Mn, N^)- 
The matrix La'' is then given by 

L2p-i'^ = iV;, L2/ = M;. (p = l,---,n) (B.7) 
In matrix notation this can be written as 

L = (iVi,Mi,---,iV„,M„). (B.8) 
To show I det L| = 1 we note that any vector ^ G Z^" can be expanded uniquely as 

i = Y,W'N, + VM^) with ap = -0(Mp,O, &p = -0(e,iVp), (B.9) 

where the coefficients and hp are all integers. In particular we consider expansion of unit 
vector s e a with (ca)^ = 5^ as 

e, = Y^^AapNp + 5,pMp), (B. 10) 



where Aap and -Bap are all integers. In terms of components we have 



(B.ll) 



These can be expressed in a matrix form as 



A 



11 



.A 



2nl 



^in Bn 

^2nn B2nl 



Bin \ 
B^nn/ 



/iVi^ 



-''n 

Mi^ 



Ari2"\ 
Mi^»^ 



Taking the determinant of both sides of this expression, we find that 

I det L\ = I det(Ari, ■ ■ ■ , Ar„, Mi, ■ ■ ■ , M„)| = 1. 

This completes the proof of the lemma. 

From ()B.1|) and ()B.13|) we obtain det u = det (j). This immediately gives 



1^1- --l^n 



\^axbx---anbn'Paxb-i ' ' ' <Pa„6„ | • 



(B.12) 



(B.13) 



(B.14) 
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c Existence of a Standard Basis 



Using the the lemma given in Appendix we can always find 2n sets of integral vectors 
Np = {Np") and Mp = {Mp") satisfying 

F{Mp, N^) = Vpbp,, F{Mp, M,) = F{Np, N^) = O, 
|det(iVi,Mi,---,iV„,M„)| = l, (Cl) 

where z/i, ■ ■ ■ , i/„ are positive integers with the property that Up divides Ug for p < q. To 
make the paper self-contained we give a proof of the existence of a set of orthonormal basis 
vectors [C^, - ■ ■ , C^") that converts the field strengths into a Standard form ()5.1|) and satisfies 
the conditions (jEUDI)- We basically follow Ref. pP. 

We first note that an arbitrary vector V G R^" can be expanded in three ways as 

V = Y, v^ea = Y, u^L, = Y, w\\ (C.2) 

a a a 

where the bases Ca-, La = (A^i, Mi, • • • , Nn.M^) and in R^"- are defined by 

(e.)^ = 6l {Laf = U\ iCf = Cl (C.3) 

This gives 

= Y,u'^L,\t- (C.4) 

a,b 

The transformation matrix 

{LR),' = YLaXl (C.5) 

c 

relates B = (Bab) with z/ = (uab) as 

u = LFL^ = LRB{LKf. (C.6) 

We then define n dimensional vector subspaces spanned by {A^i, ■ ■ ■ , A^^} and {Mi, ■ ■ ■ , M„} 
as 

M- = {V e R^'''\V = ^t^^iVp, e R}, (C.7) 

p 

M+ = {V e R^'^IV = ^t;^Mp, yP e R}. (C.8) 

p 

The field strength 2-form ()2.Hj) is degenerated on and in the sense 

F{Np, Ng) = F{Mp.Mg) = 0. (C.9) 
25 



Furthermore, M"*" and M are maximal degenerate subspaces of R^". This implies F(y, W) 7^ 
O for some K G M+ and W' ^ M+. 

The R^" can also be decomposed into n two- dimensional orthogonal subspaces Wi, ■ • • , 
Wni where Wk is spanned by the orthonormal basis vectors C^*^""^ = {Ca^~^) ^^^d C,'^^ = (Ca'^) 
satisfying 

E FM!'-' = -BlC-\ C = ^ E ^'^^Cf (C.IO) 

b,c b 

We note that the vector subspaces of Wk defined by 

= Wkn M", w+ = Wkn (m-)^ (c.ii) 

must be one-dimensional. This can be seen as follows; If dini WjT = 2, we have Wj7 = Wk C 
M~ . Since F is degenerated on M~ , we obtain F{V, V') = O for arbitrary V, V' G Wj7 . This 
contradicts the fact that the basis vectors (^^fc-i g^j^^ f^2k ^^gg^ed above satisfy 

Fie'-\n = -§^^o. (C.12) 

On the other hand if dimVF^ = 2, we see that W^ = Wk C (M")-*- is orthogonal to M~ . 
Then we can expand Np as 

Np=J2iclC-' + dlC)- (C.13) 

We now choose an arbitrary vector V G Wj^ and consider the n + 1 dimensional subspace 
M" D M" of R^" spanned by V and {A^i, ■ ■ ■ , A^„}. Using ()C.13|) . we immediately find 
F{V,Np) = O for any p. This implies that F is degenerated on M~ . This contradicts the 
fact that M~ is a maximal degenerate subspace. We thus conclude that dimW^ = 1. 
We can alway find an orthonormal basis (C^'^"^? C^'^) to satisfy 

Since PF^ is orthogonal to M~ we obtain 

5^iV/Cf-' = 0. (C.15) 

a 

This completes the proof of the existence of satisfying ()6.10|) . 
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D Computation of the Normalization Constant 



In this appendix we compute the normalization constant of the wave function ()6.45|) . Due 
to the orthogonahty relations ()6.48|) we have only to compute the integral 



D 



"'O ir, ... =-„ ^' ... ^' 



C (f^^ ^ ^ g-'rEp,,(9'"'W)(!^p€^-"p)(t',5''-"<,)+iTEp?p"p 

"'O ST, ... TT„ — / ... — / 



"^0 u. ... ^ 

d^n ^-27, Y:p,q 9"" {^pi''-rp+^ep){uqii-rq+^eq) 
■00 J —00 



k\,- - ,k„=—oo 



2 2 z/l ■ ■ ■ i/„ 



(D.l) 



This justifies the overall normalization of ()6.45|) . In two dimensions we have gu = v. The 
normalization constant appearing in ()H.17|1 is consistent with (jD.ljl . 
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